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The range of a high temperature rotating cylinder viscometer has been extended so that 
viscosity measurements can be made between 10° to 10 10 poises to within 2-percent accuracy. 
This involves three different techniques. After a calibration constant has been determined 
for the apparatus in the standard oil range, the other two constants for the higher viscosities 
were self-determined. In making measurements at the low viscosities, 10° to 10 5 poises, the 
outer cylinder is rotated at constant speed and the torque measured on the inner cylinder; 
from 10 4 - 5 to 10 7 - 5 poises the inner cylinder is rotated through an angle and timed as it returns 
to its zero position, and from 10 5 - 5 to 10 10 poises the inner cylinder is driven through an angle 
at constant torque and timed as it traverses this angle. The Fulcher equation viscosity 
values obtained from previous measurements on Standard Glasses No. 710 and No. 711 by the 
fiber elongation and rotating cylinder methods have been compared with the present results. 



1. Introduction 

In 1924, English [1] ! reviewed methods for 
measuring viscosity of glasses and concluded that 
the Margules [2] rotating cylinder technique was 
readily adaptable to this work and gave much more 
reliable and repeatable results than the previous 
techniques. His apparatus consisted of a fixed 
cup and an inner rotating cylinder driven by the 
torque produced from a counter-weight and pulley 
system. With this equipment English was able to 
make A r iscosity measurements from 10 1 to 10 7 
poises 2 on soda-lime-silica glasses. Later Parks 
et al., [3,4] using the same technique, made viscosity 
measurements on B 2 3 glass from 10 5 to 10 11 P. 

About this time Field [5] and Lillie [6], working 
with molten glasses, improved the rotating cylinder 
viscometer by eliminating the disadvantages in- 
herent in the counter-weight and pulley system. 
They rotated the cup at a constant angular velocity 
and measured the force on the inner cylinder by 
means of the angle of deflection of a torsion wire. 
This not only allowed Lillie to measure absolute 
viscosities but also gave him a higher degree of 
accuracy. This method was useful in the range 
10 to 10 5 P. 

Lillie went still further and extended the range 
of his viscometer by developing the aperiodic method 
of measuring viscosities. In this method the inner 
cylinder was deflected and allowed to decay to 
equilibrium as the outer cylinder remained at rest. 
This extended the range of his viscometer up to 
10 7 P. 

In later years, other glass investigators such as 
Babcock [7], and Robinson and Peterson [8], have 
used equipment similar to LilhVs because of the 
comparatively high accuracy. A further improve- 
ment to Lillie's apparatus was made by Bockris 
and Lowe [9], who inserted in the suspension system 



1 Figures in brackets indicate the literature references at the end of this paper. 

2 One poise is 1 g/cin sec which is equal to 0.1 kg/in sec. 



of the inner cylinder a coil and magnet for electro- 
magnetically measuring the torque on the inner 
cylinder. Their apparatus was used to make vis- 
cosity measurements on molten silicates in the range 
1CT 1 - 3 to 10+ 4 - 2 P and up to 1800 °C. 

The rotating cylinder apparatus has become, over 
the years, a standard method for measuring vis- 
cosities of glasses at high temperatures, but other 
techniques had to be developed to determine vis- 
cosities above 10 7 P. Above 10 9 P the fiber elonga- 
tion method of measuring viscosities has proven to 
be both reliable and accurate. In the gap between 
these two methods there has been relatively little 
work. 

Hampton [10] measured viscosities from 10 8 to 10 13 
P observing the rate of deflection of a beam held in 
cantilever fashion. Recently, Hagy [11] developed a 
method of determining the viscosity in the range 10 8 
to 10 15 P by the viscous bending of a simple beam 
under load. Kelley [12] et al., have developed a 
penetrometer to measure viscosity of glasses from 
lO 5,7 to 10 9 P by measuring the rate of penetration 
of a solid cylinder under load into the glass. 

Even though measurements can now be made con- 
tinuously over the whole viscosity range of glasses by 
various techniques, it seems desirable to make 
measurements with the rotation viscometer over as 
wide a range as possible since it has proven to be the 
most reliable and accurate apparatus that has been 
used in this work. 

2. Apparatus 

The geometry of this apparatus, see figure 1, is 
similar to that of Lillie [6], Bockris and Lowe [9]. 
This equipment has been previously described in more 
detail by Napolitano and Hawkins [13]. 

The outer cylinder consisted of a platinum crucible 
5 cm in diam and 10 cm high. The bob was 1.25 cm 
in diam 3.75 cm long with 45° tapered ends and was 
suspended by a hollow platinum spindle extending 
into the center of the bob to accommodate a thermo- 
couple for measuring the temperature of the glass 
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Figure 1. Wide-range high temperature rotating cylinder 
viscometer. 

before and after viscosity measurements. The cruci- 
ble was mounted on a mullite pedestal which extended 
out through the bottom of the furnace and was at- 
tached to the rotating mechanism. The rotating 
mechanism provided a range of crucible speeds from 
0.03029 rpm to 1.3493 rpm. 

An electromagnetic torque balance [9] based on the 
principle of a large moving coil galvanometer was 
used to measure the torque by a null detection 
method. The galvanometer was made up of a 
permanent magnet with a 1650 G gap field strength 
and a coil which was an integral part of the suspension 
apparatus. The coil consisted of a rectangular alu- 
minum frame with about 25 m of No. 34 insulated 
copper wire, wound around it. The resistance of the 
coil was 222 12. In series with the coil were several 
resistors which allowed continuous adjustment be- 
tween 2 kft and 5 Mil and four L and N Standard 
Fixed Resistors 10, 25, 100, and 1000 ft. The source 
of power consisted of 12-2 V cells. The balance was 
observed through a telescope on a 1 m radius scale 
located 1 m distant from a reflecting mirror set on the 
axis of the suspension shaft just below the coil. The 
amount of current through the coil was found by 
measuring the voltage across a standard resistor with 
a Type K potentiometer. Polarity in the coil would 
be reversed to accommodate reversal of rotation. 

3. Operation and Calibration 

The extension of the viscosity range of this 
apparatus was accomplished by operating it in three 
different modes. A careful examination of its 



equation of motion is very important before de- 
scribing each mode of operation. The general 
equation of motion for the viscometer is 



i^+w 



dd 
di 



(1) 



/ — moment of inertia; 

6 — angular displacement; 

t — time; 

rj — dynamic viscosity; 

Ki, K 2 , and Kz are constants of the apparatus; and 

i — current through coil. 

Thus the first term of the equation represents the 
resistance to acceleration due to inertia. The sec- 
ond term is the viscous drag. The third term rep- 
resents an elastic torque due to an angular deflection 
of the suspension, and the last term represents the 
magnetic torque applied to the system. The limi- 
tations of eq (1) are discussed in the appendix. 

Method 1. Over the lowest viscosity range 3 log 
)/=0 to log rj = 5 the outer cylinder is rotated at a 
fixed angular velocity. The viscous drag of the 
liquid produced a torque on the inner cylinder. A 
measure of this torque is made by applying an 
appropriate current to the coil such that the re- 
sulting magnetic torque will bring the inner cylinder 
back to its rest position at 0=0. 

Therefore the first term of the eq (1) is zero, be- 
cause there is no angular acceleration in the system 
when measurements are being made. The second 
term is present and significant. It is retained in the 
equation of motion for Mode I, but since the inner 
cylinder is held at rest, it is the angular velocity of 
the outer cylinder that is used in this term. The 
third term is zero when the system is balanced, since 
the angle across the suspension is zero. Even so, 
it controls the sensitivity. By changing suspensions 
the value of K 2 can be changed and with it the sensi- 
tivity. Since the fourth term is significant, it is kept. 
Thus the equation of motion for Mode I is 



or 



vKtf^Ksi 



rj=Cii/cx) 



(2) 
(3) 



where Ci=W^ and o)=-j- (angular velocity). 



To obtain G x in eq (3) the viscometer was calibrated 
at 25 and 40 °C with NBS Standard oils. These oils 
ranged in viscosity from 1 to 500 P. Some of Dow 
Coming's silicons have also proven to be a valuable 
check on the calibration constant at the higher 
viscosities. Thermal expansion corrections on C u 
known for room temperatures, have to be made for 
accurate viscosity measurements at temperatures, T, 



s During calibrations (at room temperature) measurements were made for 
viscosities below IP. At elevated temperatures the lower viscosity limit was 
10 2 P because of the characteristics of the NBS Standard Glasses and the 
temperature limit of the furnace. 
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up to 1400 °C according" to the following equation 



C 1 (T)=O l250c (l+cc[T-25])- 3 



(4) 



where a is the linear expansion coefficient of platinum. 

Table 1. Typical scries of viscosity measurements made on 
Standard Glass No. 711 by Method I 

Standard resistor 100f2; 0.030 cm diam. suspension 



Rotation 


i 


CO 


A r 


Clock wise___ __ ... 


149. 706 
149. 605 

87. 026 
86. 655 

49. 991 
49. 905 

36. 892 
36. 630 

149. 010 
148. 505 


rpm 
1. 3493 
1. 3493 

.7845 
.7845 

. 4513 
.4513 

.3332 
. 3332 

1.3493 
1. 3493 


fiA/rpm 


Counter-clockwise 

Clockwise 

Counter-clockwise 

Clockwise.. _ .__ __ 


110.91 
110.70 


Counter clockwise 

Clockwise- 


110. 68 


Counter-clockwise 

Clockwise 

Counter-clockwise 


110.34 
110.25 


Average. ______ ._. 


110.58 







Temperature before start of measurements was 1050.7 °C and after measure- 
ments was 1051.3 °C: the average being 1051.0 °C. 



A typical series of viscosity measurements made 
by this mode of operation is shown in table 1. In 
making these measurements four different speeds are 
used with readings taken in the clockwise and 
counterclockwise rotation for each speed. The 
first speed selected is always repeated at the end of 
the series to determine whether the glass viscosity 
and, therefore, the temperature, has remained essen- 
tially constant. Note the slight change in N which 
reflects the increase in temperature. From these 
measurements the viscosity is 

r ? =A^XC y 1(1050Oc) = 110.6Xl0- 6 XH.99X10 +6 = 1326P 



or 



log 17=3.123 at 7 = 1051 °C. 



Method 2. For viscosities from about log ^=4.5 
to log 77=7.5 the aperiodic method developed by 
Lillie is used. In this decay mode of operation the 
outer cylinder remains at rest while the inner cylinder 
is displaced by an angle 6 and allowed to return to its 
zero position under the force of the torsion wire. 
Since the system is an over-damped harmonic oscilla- 
tor one observes an exponentially decaying angular 
velocity as the inner cylinder approaches its equilib- 
rium position. The term involving the moment of 
inertia in eq (1) is extremely small and can be 
neglected, giving rise to the following equation : 



dd 
df 



— C 2 t 



(5) 



where 



r _K t 



Integrating 



we have 



f\ J j=-fca 

JO, V Jh 

CWr-U) 



In 61/62 



(6) 



The calibration constant (CV) for Method II can 
be calculated once the viscometer has been cali- 
brated by Method I. First an independent meas- 
urement of the current required for a given angular 
deflection of the inner cylinder is made. This gives 
us: 



K s , 



K, 



i=Cii. 



(7) 



If the current, i, is given in amperes, and the 
angle, 0, in units of 17' of arc which correspond to a 
deflection of 1 cm in a circular scale 1 m distant 
from the central axis of the spindle; one has 



Cf- 



25 



5.094X10" 



=4907.7. 



From eq (3), (5), and (7) one may obtain C 2 from 
Ci and C 4 , according to the following: 



2=JT U 



(8) 



where u changes units from revolutions per minute to 
centimeters per second. The rates in centimeters per 
second were read on the circular scale. 



u= 



cm 



min 



l x ^X100X2 =2()>95> 



Thus 



sec revolution 60 



n 20.95X1.212X10 6 _ __,. 
Ci= 4907.7 ~ 5 ' 174 - 



Table 2. Typical series of viscosity measurements made on 
Standard Glass No. 711 by Method II 

0.063 cm diam suspension 



Headings 


01 


02 


t 2 -h 


1. . 


cm 
25 
25 
25 
25 


cm 
2.5 
2.5 
2.5 
2.5 


sec 
557. 31 


2__. _ 


545.40 


3 


548.49 


Average.. ._ _ __ 


550. 40 







Temperature before start of measurements was 704.2 °C and after measure- 
ments was 704.4 °C, the average being 704.3 °C. 

A typical series of viscosity measurements made 
by this mode of operation is shown in table 2. A 
phosphor bronze suspension wire of 0.063 cm diam 
was used. Before starting, the zero was checked 
and properly adjusted. The current was applied to 
the coil so that the inner cylinder was rotated from 



451 



its equilibrium position to an angle corresponding* to 
about 30 cm on the telescope scale. The current 
was turned off and time was measured for the dis- 
placement to decay from 25 to 2.5 cm. Thus the 
viscosity is 

i.= C, ' f - ( ^ fl) = 517 ^5 ) - 4 =1.236X10'P 



or 



In 0i/0 2 2.303 

log 77=6.092 at T= 704.3 °C. 



Method 3. For viscosities from about log 77 = 5.5 
to log 77=10 this method was used. The outer cyl- 
inder again remains at rest, while a torque is applied 
to the inner cylinder by means of the magnetic field 
interaction, and the angular velocity is measured 
about 0=0. Thus the term K 2 6 is negligible and 
can be omitted. Since measurements are made in a 
steady state the acceleration is zero, and we have the 
following: 



dd K 3 i n . 



(9) 



integrating 



we have 



^2 I ^2 

y\d6= 

i Jh 



CMt 



n=C z i 



t 2 — 1\ 



(10) 



ft and d are equal, but expressed in different units: 

(7 3 =^Ci=20.95X1.212X10 6 =2.539X10 7 , (11) 

where u is the same conversion factor used in cal- 
culating C 2 . 



Table 3. Typical series of viscosity measurements made on 
Standard Glass No. 710 by Method III 



Standard resistor 10ft 


0.063 cm 


diam suspension 


Rotation 


02-01 


h-U 


% 

AXKH 


Clockwise.. _ _ _. 


cm 
1 
1 
1 
1 
1 


sec 
197. 58 
195. 48 

193. 62 
191. 42 

194. 52 


196. 13 


Counter clockwise 

Clockwise-. 


196. 13 
196. 13 


Counter clockwise 

Average— 


196. 13 
196. 13 







Temperature before start of measurements was 706.0 °C and after measurements 
was 707.6 °C: the average being 706.8. 



A typical series of viscosity measurements made by 
this mode of operation is shown in table 3. A 
measured amount of current is applied to the coil, 
which rotates the inner cylinder at a constant angular 
velocity. The time required to traverse 1 cm about 
the equilibrium position (6=0) on the scale is 
measured. By appropriately selecting the current, 
measurements can be made up to log 77=10 P within 
reasonable waiting periods. 



From these measurements the viscosity is 

C^h-U). 



-i=2539Xl0 4 X 194.52 



02-0i 

K10X196.13X10" 4 

=9.6865X10 7 or log 77 = 7.986 at T= 706.8 °C. 

4. Discussion of Results 

Measurements were made on NBS Standard 
Viscosity Glasses Nos. 710 and 711. These results 
are given in tables 4 and 5. Only the measurements 
made by Methods II and III are listed here, since 
Method I and the fiber elongation methods of 
measuring viscosities have been covered in previous 
work [13, 14]. 



Table 4. Viscosity values of Standard Glass No. 710 by 
Methods II and III compared with values calculated from 
Method I and fiber elongation method* 



Temp. °C 


log ?7ob a 


lOg 77calc 

Fulcher 




Method II 


Method III 


equation 


951.7 
902.3 
851.2 
810.3 
774.7 
773.4 
757.5 
757.1 
725.9 
706.8 
693.2 
677.8 
657.9 
636.0 


4.581 

5.028 
5. 629 
6.169 

~~6~999~ 


~~6.~718~ 
6.734 

_ ~7~6l3~ 

7.608 
7.986 
8.298 
8.675 
9.178 
9.827 


4.553 

5.033 
5.614 
6.158 
6.702 
6.723 
6.993 
7.000 

7. 585 
7.984 

8. 290 
8.660 
9.182 
9.821 



*Fulcher equation constants taken from Napolitano and Hawkins (Lab. A) [13 

Table 5. Viscosity values of Standard Glass No. 711 by 
Methods II and III compared with values calculated from 
Method I and fiber elongation method* 



Temp. °C 


lOg rjoba 


lOg 77calc 

Fulcher 




Method II 


Method III 


equation 


774.2 
770.5 
755.0 
734.4 
732.0 
725.0 
704.3 
677.6 
651.7 
631.5 
607.4 
574.7 
554.9 
528.1 


5.220 
5.255 
5.438 
5.696 
5.722 
5.831 
6.092 
6.488 


"~5.~436" 

5.690 
5.713 
5.813 
6.082 
6.480 
6.909 
7.281 
7.752 
8. 423 
8.904 
9.653 


5.229 
5.270 
5.446 
5.695 
5.725 
5.815 
6. 092 
6. 482 
6.900 
7.256 
7.723 
8.454 
8.924 
9.674 



*Fulcher equation constants from Napolitano and Hawkins (Lab. A) [14]. 

The viscosity values obtained by the rotation 
method (Method I) and the fiber elongation method 
were fitted by a least square calculation to the Fulcher 
[15] equation. Fulcher stated that an empirical 
equation of the form 



log rj=A+B!(T- To) 



(12) 
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accurately relates viscosity to temperature in many 
glasses. These equations become: 

log 77= -1.626 +423 9/(r- 265.7) 

forSt'd. Glass No. 710; (13) 
log T7=-1.654+4317/(r-- 147.0) 

forSt'd Glass No. 711. (14) 

These are shown in figure 2 as the line curves in 
order to illustrate the viscosity values obtained in 
this study, and to compare the results from Methods 

II and III with previous measurements. The vis- 
cosity values obtained by these methods for the two 
standard glasses fall so close to the curves obtained 
from the previous measurements, that no deviation 
can be seen. In order to inspect these deviations 
more closely, the observed viscosity values were 
compared with the values calculated from eqs (13) and 
(14), for each temperature. These deviations are 
plotted in figure 3 as well as those from Method I 
and from the fiber elongation method. The gap in 
rotation and fiber elongation measurements has been 
bridged. Using all four techniques one can make 
measurements over a continuous viscosity range from 
10 1 to 10 12 P. The scatter of the viscosity values 
obtained by the driving technique is less than that 
obtained from the fiber elongation technique. 

In figure 3, the zero lines represent the Fulcher 
equation as calculated from Methods I and IV. Thus 
the agreement between the data from Modes II and 

III and this line is a test not only of the precision but 
also of the accuracy of these measurements. Since 
the errors fluctuate to both sides of the zero line, any 
systematic deviation should be small. In fact the 
systematic deviations were calculated for each mode 
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IIOO IOOO 900 
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TEMPERATURE-°C 
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Figure 2. Plot of log viscosity versus temperature of NBS 
Standard Glasses No. 710 and 711. 

Circles represent data taken by driving inner cylinder while triangles represent 
data taken with the decay mode of operation. The line curves represent calcu- 
lated values from Napolitano and Hawkins (Lab. A) [h.m]. 



using all the measurements taken in both standard 
glasses. Their values were 0.30 percent for Mode II 
and 0.45 percent for Mode III. The standard devia- 
tions were calculated similarly. They were 0.010 log 
viscosity (2.3%) for Mode II, and 0.013 log viscosity 
(3.1%) for Method III. 

Thus the equation of motion of the apparatus 
(eq (1)) has been proven correct as well as the calcu- 
lated relations between the various calibration con- 
stants. This has been indicated in the appendix. 
The relatively little loss in accuracy due to the exten- 
sion of the viscometer is probably due to the larger 
temperature dependence of the viscosity at lower 
temperatures. 

5. Summary 

The range of the rotating cylinder high temperature 
viscometer was extended up to 10 10 P without loss of 
accuracy. Measurement at the higher viscosities 
could be made within reasonable time periods. 

This extension was accomplished by operating the 
viscometer in three different modes: (1), rotating the 
outer cylinder at uniform angular velocity and meas- 
uring the torque in the inner cylinder; (2), the outer 
cylinder remains at rest while the inner cylinder is 
displaced by an angle and allowed to return to its 
zero position under the force of a torsion wire; and 
(3), the outer cylinder again remains at rest while a 
torque is applied to the inner cylinder by means of the 
magnetic field interaction and the angular velocity is 
measured about 0=0. The viscometer was cali- 
brated by Method I, and from this the constants for 
the other two methods were calculated. This cali- 
bration was checked with NBS Standard Samples 
Nos. 710 and 711, and found to be within 0.5 percent, 
well inside the statistical fluctuation, whose standard 
deviation was 0.010 log viscosity (2%) for Methods 
I and II and 0.013 log viscosity (3%) for Method III. 
This is considerably better than that obtained for the 
fiber elongation data 0.034 log viscosity (8%). 



6. Appendix 

The second term of eq (1) assumes perfect laminar 
flow, which is not always obtainable in a rotation 
viscometer. Three other flow patterns are possible 
in such a system — Taylor vortices [16], convection 
currents and turbulence. 

Taylor has made an extensive study of nonlaminar 
flow of viscous liquids between rotating cylinders and 
some of his conclusions can be applied to the geom- 
etry of the present apparatus. Taylor found that 
his vortices appeared above a critical value of co/77. 

The largest value of — for Mode I in the present 
V 
investigation was 0.52. This compares with Taylor's 
experimental critical value much greater than 1000 
and calculated at infinite speeds. For Mode II they 
can also be neglected since its highest u is 0.1 rad/sec 

and lowest rj is equal to 5X10 4 giving— =5X10~ 5 
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Figure 3. Discrepancies between the observed log viscosity values and those calculated from the Fulcher equation for NBS Standard 

Glasses No. 710 and No. 711 are plotted versus log viscosity. 
Open triangles — fiber elongation; closed circles— driving inner cylinder; closed triangles— driving outer cylinder; and open circles— rotating outer cylinder. 



compared with Taylor's experimental value of 30. 
Similarly, they are absent for Mode III. One 
should note that the ratio of diameters of the present 
apparatus is considerably larger than that of Taylor's 
and also that the center cylinder is very much shorter. 

However, since the values of — used in this work are 

many orders of magnitude smaller than Taylor's 
critical values above which vortices appear, one is 
justified in assuming the absence of said vortices in 
this apparatus. 

Following similar reasoning one can rule out the 
effects of convection currents in molten glasses when 
the cup has a temperature gradient of less than 5 
deg. In a previous paper [13] it has been shown 
that the temperature gradients in the cup are much 
less than this. Since turbulence requires an even 

higher value of — , it definitely is not present. 
1 



The expectation that nonlaminar flow does not 
exist under the experimental conditions used was 
checked experimentally by plotting the torque versus 
angular velocity for both Modes I and III. The 
relationship was found to be linear. In Mode II a 
plot of displacement versus time was found to be 
exponential. This is direct experimental evidence 
that these nonlaminar effects are negligible, in agree- 
ment with the analysis given above. 

Two other possible sources of error exist in our 
apparatus. They are (1), any effects due to the bob 
and cup not being perfectly centered; and (2), the 
dependence of end effects on the depth of immersion 
and clearance from the bottom of the cup. 

Due to the large ratio of cup-to-bob diameters, 
precise centering is difficult to achieve. However, 
this difficulty is more than compensated for by the 
fact that slight departures from centering are not 
critical. After reasonable care is taken to center 
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the cup and bob, a quantitative but sensitive measure 
of the alinement obtained can be determined from the 
stability of the null position in making viscosity 
measurements according to Mode I. The results 
thus obtained are applicable to Modes II and III. 

Lillie [17] and Rait [18] have shown that by having 
conical ends on the bob the dependence of end effects 
due to bob clearance is considerably reduced. This 
has been confirmed by the authors, who raised and 
lowered the cup 3 mm during viscosity measurements 
with negligible effects. 

In conclusion it has been shown that eq (1) is 
valid for viscosities between 10° and 10 10 P. 
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